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ABSTRACT 

Selection of active galactic nuclei (AGN) in the infrared allows the discovery of AGN whose optical 
emission is extinguished by dust. In this paper, we use the Spitzer Space Telescope First Look Survey 
(FLS) to assess what fraction of AGN with mid-infrared luminosities comparable to quasars are missed 
in optical quasar surveys due to dust obscuration. We begin by using the Sloan Digital Sky Survey 
(SDSS) database to identify 54 quasars within the 4 deg 2 extragalactic FLS. These quasars occupy a 
distinct region in mid-infrared color space by virtue of their strong, red, continua. This has allowed 
us to define a mid-infrared color criterion for selecting AGN candidates. About 2000 FLS objects 
have colors consistent with them being AGN, but most are much fainter in the mid-infrared than 
the SDSS quasars, which typically have 8/xm flux densities, Sg.Oj ~ lmJy We have investigated 
the properties of the 43 objects with Ss.o > lmJy satisfying our AGN color selection. This sample 
should contain both unobscured quasars, and AGN which are absent from the SDSS survey due to 
extinction in the optical. After removing 16 known quasars, three probable normal quasars, and eight 
spurious or confused objects from the initial sample of 43, we are left with 16 objects which are likely 
to be obscured quasars or luminous Seyfert-2 galaxies. This suggests the numbers of obscured and 
unobscured AGN are similar in samples selected in the mid-infrared at Ss.o ~ lmJy. 
Subject headings: quasars:general - galaxies :Seyfert - infrared: galaxies 



1. INTRODUCTION 

The number of AGN missed in optical or soft X-ray 
surveys due to obscuring columns of dust and gas in the 
host is still an open question (e.g., Webster et al. 1995; 
Cutri et al. 2001; Gregg et al. 2003; Norman et al. 2002; 
Richards et al. 2003). Understanding and quantifying 
the population of obsured AGN is important if we are to 
obtain a complete picture of the build-up of black holes 
in the nuclei of galaxies. The close link between bulge 
mass and black hole mass (e.g. Magorrian et al. 1998) 
implies that understanding this history will also help us 
to understand the galaxy formation process as a whole. 

The local mass density in black holes is dominated 
by the ~ 1O 8 ' 5 M black holes in L* elliptical galaxies. 
When accreting at Eddington rates, these correspond to 
luminous quasars (Mb ~ —25). Recent work on the 
X-ray background, however, indicates that most of the 
mass build-up may occur on longer timescales at sub- 
Eddington rates, mostly in obscured AGN (Cowie et al. 
2003). Besides X-ray searches, other techniques for find- 
ing obscured AGN have also been employed. A sample 
of luminous narrow-line AGN has been selected from the 
SDSS (Zakamska et al. 2003). Quasar samples selected 
in the near-infrared from the Two Micron All-Sky Survey 
(2MASS) (Cutri et al. 2001; Glikman et al. 2004) have 
identified a number of quasars with significant optical 
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extinction and X-ray columns. Even in the mid- infrared 
(MIR) , a highly-obscured AGN will still have significant 
extinction. However, the disk or torus-like geometry of 
the obscuring material means that even if the extinction 
to the broad-line region is high, MIR light will usually 
be visible, unless the geometry is truly edge-on, or the 
AGN is completely enveloped in a dense cloud of dust and 
gas (Pier & Krolik 1993; Efstathiou & Rowan-Robinson 
1995). Thus an MIR survey can usefully address the 
question of the number of obscured AGN. 

Studies with the Infrared Space Observatory (ISO) 
have shown that the strong MIR continuum associ- 
ated with AGN provide a unique spectral signature that 
can be used to distinguish AGN from starbursts. The 
MIR continuum from galaxies arises mostly from three 
sources: Hn regions dominated by emission from very 
small dust grains (producing a steeply-rising continuum 
at 12-16/im), photodissociation regions dominated by 
bands of Polycyclic Aromatic Hydrocarbon (PAH) emis- 
sion, and AGN dominated by a strong 3-10/im contin- 
uum. Laurent et al. (2000) show that the spectral en- 
ergy distributions (SEDs) arising from each of these is 
sufficiently distinct to allow discrimination of AGN from 
star-forming galaxies based on their MIR SEDs. Using 
these ideas, Haas et al. (2004) combined near- and mid- 
infrared colors to select AGN from an ISO parallel sur- 
vey. With the advent of Spitzer (Werner et al. 2004) 
it has become possible to obtain MIR photometry for 
large samples of field galaxies, and thus to use the MIR 
SEDs to identify specific populations of AGN and star- 
forming galaxies. In this paper we use SDSS quasars to 
provide an empirical localization of the AGN population 
in mid-infrared color space. We then select a sample of 
candidate obscured AGN with 8/im flux densities (5*8. o) 
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> lmjy, and compare their optical identifications and 
estimated redshifts to those of the SDSS quasars with 
5*8.0 > fmJy. We use the main field of the extragalac- 
tic component of the Spitzer First Look Survey (FLS), 
a shallow, 4 deg 2 survey with the Infrared Array Cam- 
era (IRAC; Fazio et al. 2004) and the Multiband Imag- 
ing Photometer for Spitzer (MIPS; Rieke et al. 2004), for 
which an extensive multiwavelcngth ancillary dataset ex- 
ists. We assume a cosmology with 51m = 0.3, J7a = 0.7 
and H = 70kms- 1 Mpc- 1 . 

2. THE SPITZER FLS MAIN FIELD DATASET 

The FLS observations were made in December 2003 
(program ID 26, astronomical observation request (AOR) 
IDs 3861504, 3861760, 3862016, 3862016, 3862272, 
3862528, 3862784, 3863040, 3863296, 3863552). We have 
used preliminary versions of the IRAC catalog (Lacy 
et al. 2004) and the MIPS 24/zm catalog (Fadda et al. 
2004b). The IRAC catalog has flux density limits (5a in 
a 5 diameter aperture) in the four IRAC channels with 
nominal central wavelengths of 3.6/im, 4.5/zm, 5.8/im and 
8.0^m, of S 3 . 6 w 7/iJy, 6*4.5 ~ 8/iJy, 5 5 . 8 « 60/xJy, and 
Ss.o ~ 50^iJy. The MIPS 24^m catalog has a flux density 
limit of « 300/iJy. 

3. THE SDSS QUASARS 

The Sloan Data Release 1 (DR1) quasar survey 
(Schneider et al. 2003) contains 54 quasars which fall 
within the main field of the FLS. All are detected with 
IRAC, and all but one at 24^im with MIPS. Three more 
quasars were found as a result of spectroscopic follow-up 
of the radio survey of Condon et al. (2003) using the Kast 
Spectrograph at Lick Observatory. 

4. THE POSITION OF SDSS QUASARS IN MIR 
COLOR-COLOR PLOTS 

An IRAC color-color plot, using all four broad-band 
channels of the IRAC instrument, is shown in Figure 
1. The dots indicate the location of w 16000 objects in 
the main field catalog. Plotting the 8.0^m/4.5^m ratio 
against the 5.8/j.m/3.6/im ratio makes the color-color plot 
effective at separating objects with blue continua from 
those with red. Most objects have blue colors in both 
axes. These are most likely stars, or low-redshift galaxies 
whose SEDs are weak in non-stellar light, e.g., elliptical 
galaxies. From this clump two distinct "sequences" can 
be seen. One has blue colors in 5*5.8/5*3.6 and very red 
colors in Sg.o/S^. These are most likely low-redshift 
(z < 0.2) galaxies with the centers of their 6.2 and 7.7/xm 
PAH emission bands (Puget & Leger 1989) rcdshiftcd 
into the IRAC 8.0£im filter. The other sequence has red 
colors in both pairs of filters and it is on this sequence 
that the SDSS and radio-selected quasars lie. Note that 
strong stellar light will shift the 55.s/53.6 ratio to the 
blue, as illustrated by the two filled squares in Figure 
1. These objects are classed as quasars in the SDSS 
database, but have resolved host galaxies in the IRAC 
images, and so are classed as Seyfert-ls in this study. 

5. MIR SELECTION OF CANDIDATE AGN 

We used Spitzer colors of the SDSS quasars to empiri- 
cally define a color selection which will pick out objects 
with MIR SEDs consistent with their being AGN. This is 
the region enclosed by the dashed line in Figure 1. Our 
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Fig. 1. — An IRAC color-color plot using the main field FLS 
data. Dots represent the f» 16000 objects with "clean" detections 
in all four IRAC bands. Crosses indicate the colors of all the SDSS 
and radio-selected quasars, squares the SDSS Seyfert 1 galaxies 
and triangles the bright (Ss^m > lmjy) sample of obscured AGN. 
The dashed line shows the color criteria used to pick out the AGN 
sample. 
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Fig. 2. — An IRAC-MIPS 24/^m color-color plot. Dots represent 
objects with "clean" detections in all four IRAC bands matched to 
the MIPS 24^im catalog of Fadda et al. (2004b), totalling « 6000 
objects. Crosses indicate the colors of SDSS quasars, stars the 
SDSS Seyfert 1 galaxies and triangles the Sg.o > lmjy sample of 
obscured AGN candidates. 



modelling of MIR SEDs based on ISO spectra (Sajina, 
Lacy & Scott 2004) suggests that most of the objects 
in this region are indeed dominated by AGN emission, 
though we expect some contamination by star-forming 
galaxies close to the boundaries. There are « 2000 
sources in this region, which are therefore likely to con- 
tain AGN. Most of them are fainter in the MIR than the 
SDSS quasars, which are typically much brighter than 
the FLS flux density limits. 

To pick a sample of AGN selected in the MIR which we 
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could compare directly to the SDSS quasar sample, we 
examined the distribution of Ss.o for the SDSS sample. 
(8.0^im is the longest IRAC wavelength, and therefore 
least affected by dust obscuration.) This distribution 
showed a peak at Ss.o ~ ImJy, and we therefore as- 
sumed that at Ss.o > ImJy most of our MIR-selected 
AGN would appear in the SDSS sample if they were not 
obscured by dust. 

There are 43 objects within the AGN selection region of 
the color-color plot with Ss.o > ImJy and "clean" detec- 
tions in the IRAC FLS catalogue [i.e. unblended (about 
80% of the 4-band detections), and outside the haloes 
of bright stars and unaffected by multiplexer bleed or 
pulldown (97.6% of the survey area)]. The steps used to 
determine which of these were candidate obscured AGN 
were as follows: (1) we matched the sample of 43 bright 
AGN candidates to the SDSS quasar list, finding that 14 
of our 43 objects are known SDSS quasars; (2) we next 
matched to our list of radio-selected quasars, finding two 
objects in common; (3) we examined the objects by eye 
on the FLS mosaics, and rejected a further seven objects 
from the sample on the basis of them being either satu- 
rated stars, or confused in the IRAC images. This left 
20 objects in the sample, which are our bright, obscured 
AGN candidates. These are shown as triangles in Figure 
1, and listed in Table 1. As a check on the AGN nature 
of these candidates, in Figure 2 we plot a color-color plot 
which extends the wavelength range to MIPS 24/im. The 
large wavelength difference between 8 and 24/xm makes 
the interpretation of the MIPS/IRAC color-color plots 
more complicated, but it can be seen that most of the ob- 
scured AGN candidates lie in the same region of the plot 
as the SDSS quasars, consistent with them having SEDs 
dominated by emission from hot dust around the AGN. 
Starburst galaxies typically have much redder S24/S5.S 
colors due to a lack of high temperature dust emission 
from the AGN. 

All our 20 bright, obscured AGN candidates (and all 
the SDSS quasars) were identified on the i?-band images 
of Fadda et al. (2004a) (limiting R w 25.5). These iden- 
tifications allowed us to further refine our selection. Out 
of the 20 obscured AGN candidates, we found that three 
had the MIR to optical color and stellar morphology of 
normal quasars (one apparently behind a spiral disk), 
one candidate is very bright (R = 15.7), so is probably 
either a peculiar star, or a foreground star within ~ 1 
of a background AGN. Of the remaining 16 objects, 14 
are extended in i?-band, and thus most likely galaxies 
whose AGN is completely obscured in the optical. Two 
are faint in i?, but point-like and so are probably only 
partly obscured AGN. Figure 3 shows a plot of the ra- 
tio of i?-band flux density (So. 65) to Ss.o plotted against 
lg(Ss.o) for all 20 objects. 

6. PHOTOMETRIC REDSHIFTS OF THE CANDIDATE AGN 
POPULATION 

Of the 14 extended (galaxy) detections in i?-band, 12 
are detected on the SDSS DR1 images (which are sig- 
nificantly deeper than the SDSS spectroscopic quasar 
and galaxy surveys), and we have used their ugri magni- 
tudes to make a rough photometric redshift estimates us- 
ing Hyperz (Bolzonella, Miralles & Pello 2000), with the 
Bruzual & Chariot (1993) templates. (We exclude the 
z-band due to the likelihood of contamination by AGN 



TABLE 1 

Candidate AGN not in SDSS DR1 with Ss.o > ImJy 
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densities at 1.4GHz in the survey of Condon et al. (2003); (3) 
photometric rcdshifts, z p hoti those in brackets are obtained by as- 
suming a galaxy luminosity of L* , others are from Hyperz using 
the SDSS ugri magnitudes; (4) optical classification in the ij-band 
image; (5) this is a stellar object within the isophotes of a spiral 
disk, we assume it is a background quasar which happens to be 
behind a disk galaxy; (6) this object may be a star, see text. 





1 1 1 1 1 1 1 1 

® 


i | i i 
ffl 






H 

ffl H 








i ® 

ffl 

ffl ffl 

X 

a e x 


ffl 

a 

X 

X 


X 




>k » x x 
a x 

X 




s> 







3 3.2 3.4 3.6 3.8 

ig(s a „/VJy) 



Fig. 3. — Ss.o to -R-band flux ratio as a function of Ss.o- SDSS 
quasars are shown as crosses, point-like identifications of AGN can- 
didates are shown as triangles within circles, and triangles within 
squares denote galaxy identifications. 



emission, but believe the other bands should be largely 
free of significant contamination given the extended ap- 
pearance of the galaxies in the i?-band data) . The galax- 
ies have a mean luminosity of ~ L* , and inspection of the 
i?-band image shows they appear to be mostly early-type 
galaxies. Our photometric redshifts are only approxi- 
mate, but the magnitudes of the galaxies are such that 
few are likely to be at z > 0.7. 
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7. DISCUSSION 

Although wc emphasize that the candidate obscured 
AGN we have identified from their MIR colors require 
spectroscopic confirmation as AGN, our study suggests 
that about 50% of AGN selected in the MIR at these flux 
levels may be sufficiently obscured by dust in the optical 
for them to be missed from the SDSS quasar survey. This 
excludes AGN which are not energetically dominant in 
their host galaxies (which would plot elsewhere on the 
IRAC color-color diagram), and any AGN which remain 
obscured in the MIR. 

The median photometric redshift of our candidate ob- 
scured AGN is w 0.36 (assuming the faint identifications 
with stellar morphologies to lie at high redshift), com- 
pared to w 0.69 for the SDSS quasars with S$.o > lmJy. 
There is thus a hint that the obscured AGN may have 
lower redshifts and luminosities. Selection effects are un- 
likely to account for this difference. MIR selection en- 
sures that the AGN are selected independent of their op- 
tical properties. Whilst we might be biased towards ob- 
jects with more dust close to the AGN, this would not in- 
troduce a redshift bias if the fraction of such objects was 
constant with both luminosity and redshift. A lower me- 
dian luminosity for the obscured AGN would, however, 
be consistent with the "receding torus" model (Lawrence 
1991) in which a higher fraction of low-luminosity AGN 
are obscured, and might also be consistent with some 
studies of the X-ray background which find large num- 
bers of low luminosity, relatively low redshift obscured 
AGN (Gilli et al. 2001; Gandhi & Fabian 2003; Steffen 
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et al. 2003). 

If unobscured, our AGN would be moderately lumi- 
nous, with B-band absolute magnitudes, Mb, ~ —22, 
on the boundary between Seyfert-1 galaxies and quasars. 
The brightest, z ~ 0.6, objects would have Mb ~ —23, in 
the quasar regime, and close to the median luminosity of 
the SDSS sample. Pursuing this work to lower 8^m flux 
densities, and using spectroscopy to confirm the AGN 
nature and redshift distribution of our objects, will help 
us link the infrared to the X-ray studies. This will allow 
us to derive a consistent story for the nature and evolu- 
tion of obscured AGN, and thus for how and when the 
supermassive black holes in galaxies today accreted most 
of their mass. 
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